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Abstract 

By means of the asynchronous cellular automata algorithm we study stationary 
states and spatial patterning in an SIS model, in which the individuals’ are 
attached to the vertices of a graph and their mobility is mimicked by varying 
the neighbourhood size q. The versions with fixed q and those taken at random 
at each step and for each individual are studied. Numerical data on the local 
behaviour of the model are mapped onto the solution of its zero dimensional 
version, corresponding to the limit q -A +oo and equivalent to the logistic 
growth model. This allows for deducing an explicit form of the dependence of 
the fraction of infected individuals on the curing rate 7 . A detailed analysis of 
the appearance of spatial patterns of infected individuals in the stationary state 
is performed. 

Keywords: epidemiology, cellular automata 
2000 MSC: : 92D30, 37B15, 92C60 


1. Introduction 

As the spread of a severe disease affects large number of individuals its deep 
understanding calls for the use of mathematical models of complex dynamical 
systems. In the models of the Kermack-McKendrick type [ij, the population is 
split into characteristic groups (compartments) such as: susceptible to infection; 
infected; removed (recovered). The corresponding fractions are traditionally de¬ 
noted as S, I and R, respectively. Their evolution is obtained from a system 
of differential equations which take into account the basic mechanisms of cur¬ 
ing and infecting. Depending on the compartments involved the models are 
traditionally abbreviated as SIR and SIS (for the case of non-immune disease 
where R group is absent). Further refinements are made by introducing such 
compartments as latently infected E, or by splitting existing compartments into 
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sub-compartments based on such differentiations as age, sex, immunity level 
BSi- Usually, such models do not take into account the spatial localization 
of the individuals involved in the process, which is adequate only if the described 
populations are well-mixed and hence behave globally in space. That is why, 
such models are often referred to as zero-dimensional. 

A more advanced description should, however, be based on the local prop¬ 
erties of the system and on the deduction of the global properties from the 
local ones. This includes describing the appearance of spatial patterns (cluster¬ 
ing). Along with analytic methods spatially-structured models are also studied 
by means of the cellular automata (CA). Here one considers a system of finite 
collections of individuals. Each individual is attached to a vertex of an under¬ 
lying graph. The evolution of the system is then governed by a specific update 
algorithm which changes the state of a particular individual according to the 
states of its neighbours in the graph. This type of evolution is discrete in space 


and time and may involve explicit or implicit account for the mobility of its 
constituents sssiiEa, fni m . [isl . [m [TsI , [Til . Boccara and Cheung 
0 considered SIR CA model with two subrules: for explicit individuals mo¬ 
bility and for the infecting/curing update. They studied the influence of the 
degree of mixing due to the mobility on the epidemics spread and found, in 
particular, that for an infinity degree of mixing the time evolution turns into 
that predicted by the zero dimensional SIR model. Then the same authors 
extended their study to the case of CA SIS model [1]. A kind of phase transi¬ 
tion was found between the endemic and disease-free states, where the role of 
the order parameter is played by the stationary value of I, dependent on the 
model parameters. It was found that the mobility essentially determines the 
behaviour of the model. In particular, it reduces to that of the zero-dimensional 
SIS model in the limit of large number of tentative moves, see also Q. Another 
way of accounting mobility is as follows. The individuals are set static and their 
mobility is mimicked by varying the neighbourhood size q. This approach is 
used in the current study in two versions: (a) the neighbourhood size is fixed; 
(b) it is random at each time step and for each individual. The latter allows 
for mimicking also the spatial heterogeneity of infectivity as studied in 1^ [l8 1. 
In the study of dynamical systems, one usually finds the stationary state and 
then describes how the system approaches this state. Periodic temporal mod¬ 
ulations of parameters corresponding to some epidemic features are considered 
in Ref. 0 for CA SIS model, it was found that the inclusion of the vital 
dynamics introduces no relevant changes in its behaviour. 

In the present work, we consider CA SIS model on geometric graph obtained 
from by assigning neighbours to a given k G according to the rule: k' ^ 
k whenever \k' — k\ < for a given Rg > 0, where Rq is such that the 
neighbourhood size is equal to q. The latter paremeter can be either fixed or 
set random. The aim of the study is 


• for a fixed q, studying of the ratio of infected individuals in the stationary 
state and its dependence on the model parameters; 

• the same for globally bounded random q analysing the role of randomness; 
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Figure 1: Flow chart for the zero dimensional SIS model. 


• spatial patterning and dynamics of approaching the stationary state. 

In Sec. [21 we present the solution for the zero dimensional SIS model in the 
form used in the remainder of the paper. In Sec.jSl we consider our main model 
with the fixed neighbourhood and relate it to that of the zero dimensional version 
with rescaled infectivity. In Sec. SI we study version of the model in which the 
neighbourhood of each vertex is random of globally bounded size. In Sec. |5l 
we discuss the spatial patterning and the time evolution towards the stationary 
state. Conclusions are made in Sec. El 


2. Zero dimensional SIS model 

The flow chart for the zero dimensional SIS model [H which we consider is 
shown in Fig.[Tl Its evolution is governed by the differential equations 

S = ^I-PIS 

j = -yj + pis, ^ ^ 


where dot stands for the time derivative and /3 and 7 are infecting and curing 
rates, respectively. By the condition S = 1 — I the system turns into 


7 = (/3 - 7)/ - 




1 -7//3J ’ 


( 2 ) 


which is the logistic growth equation In the sequel, as a reference model 

we use the reduced version of m corresponding to the choice 


/3 = 1 - 7 > 7 , 


(3) 


which amounts to considering 7 restricted to 7 < 0.5 =: 7 c. In this case, 

1 - l/lc 


I{oo) = 


1-7 


and the explicit solution is given by 


m = 


/(O)J(oo) 


J( 0 ) + [/(oo)-/( 0 )]e-‘A 


(4) 


(5) 


where r = 1/(1 — y/yc) is a characteristic time for the system to reach the 
stationary state. 
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3. Cellular automaton SIS model with fixed neighbourhood size 

We consider the CA SIS model on the following geometric graph. The 
vertex set is Z^, and, for a fixed Rq > 0, the neighbourhood of a given A: S is 
defined according to the rule: k' ^ k whenever |A:' — A:| < The choice of Rg 
is made to obtain the neighbourhood size equal q. To each k there is attached a 
individual, the state of which, Sk, is either 0 (susceptible) or 1 (infected). The 
evolution of the system of N such individuals is run according to the algorithm 

• select individual k at random; 

• if Sfe = 0 do nothing; 

• if Sfe = 1 then set Sfe = 0 with probability 7 and set Sk> = 1 for one of its 
randomly chosen neighbours with probability 1 — 7; 

• perform N attempts described as above to complete a single time-step. 

It is asynchronous as far as the flip of Sk is immediate and it affects the updates of 
other individuals before the time-step is complete. On average, the probability 
p(l —> 0) for infected individual to be cured and for the susceptible individual 
to be infected, p{0 —>■ 1 ), are equal to: 

p(l ^0) = 7, (6) 

p(O^l) = {l-"f)ik/q, 


where ik is the number of infected individuals in the neighbourhood of fcth 
individual, therefore the infecting is now local. 

In the relation to this, one can refer to the following specific cases: q = 4 
corresponds to CA SIS with the von Neumann neighbourhood, referred also 
to as the contact process on a square lattice 2 l|; g = 8 realises the so-called 
Moore neighbourhood; at q = N the infecting range spans over all the system 
representing the zero dimensional SIS model © considered in the previous 
section. 

All simulations in the current study are performed on the square of 256 x 256 
individuals, total number of individuals is N — 65536. The periodic boundary 
conditions are applied in both directions. Typical times for the system to reach 
the stationary state vary considerably depending on the value of 1 ( 0 ) and on 
the proximity to the critical value 7 c, where the relaxation time is the longest. 
The value of 7 c depends on q and is defined by the condition that the fraction of 
infected individuals vanish. Near 7 c, I{oo) ^ 0.1, hence the choice of /(O) = 0.1 
reduces drastically the time needed by the system to reach the stationary state 
when 7 Ri 7c. This choice for 7(0) is used in our study and the simulations of 
2000 time-steps are performed at each 7 ranging from 7 = 0 to 7 ~ 7c. To avoid 
locking the system in the state with I{t) = 0 , we allow for at least one infected 
individual in the system. The dependence of I (00) vs curing rate 7 at various q 
is shown in Fig. [2] At q = N, it follows the solution ([5])- For q < N, the general 
shape is the same but with reduced critical curing rate 7 c. At <7 = 4, we obtain 
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Figure 2: Fraction /(oo) of infected individuals in the stationary state vs curing rate 7 for the 
CA SIS model with fixed neighbourhood size q. Symbols: simulation data; dashed curves: 
Eq. ([Sll using power law fit unii; black solid curve: Eq. (SJ for the zero dimensional SIS 
model. Values for critical curing rate 7 c are: 0.378 (q = 4), 0.408 (q = 8 ), 0.447 (g = 20), 
0.482 (g = 96), and 0.5 {q = N). 


0.368, which agree well with the value Ac = 1.64872(3) found by Sabag et 


22 |, where Ac = (1 - 7c)/7c 


7c ' 
al. 

The CA SIS model is characterised by the local infectivity ik/q ©, m 
contrary to the zero dimensional SIS model, in which case the infectivity is 
global and equal to I. Using the simulation data obtained at q < N, the zero 
dimensional SIS model can be extended to the case oi q < N hy the following 
substitution for the infecting rate: P ^ 13', where 


3' = (3 


{ik/qk) 


If the averaging is performed in the stationary state, then the dependence of 
the scaling factor for (3' on I is eliminated. Taking also into account (jS]), one 
obtains 


/ 3 '=/ 56 »( 7 ), e{-f) 


ii^k/Qk} stat 

I 


(7) 


where “stat” indicates averaging in the stationary state. By solving ([5]) with 
the substitution of /3 by P' according to 0 and taking into account © , one 
obtains I{oo) for the extended zero dimensional SIS model 


I{oo) 


1 - [1 + 61 ^( 7)]7 

1-7 


7 < 7c- 


It yields the following equation for the critical curing rate "fc- 


^(7c) 


7c 

1 -7c 


( 8 ) 

(9) 
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Figure 3: Scaling factor 6 ( 7 ) in the stationary state vs curing rate 7 for the extended zero 
dimensional SIS model with fixed q. Symbols: simulation data, dashed curves: analytic 
expression, Eq. IIIOII . Fitted values for the exponent <5 are: 0.400 {q = 4), 0.319 {q = 8 ), 0.258 
{q = 20) and 0.175 {q = 96). 


which is not equal to 0.5 now. 

To made the extended zero dimensional SIS model self-contained one needs 
certain analytic form for 0(7). It can be obtained by fitting the simulation data 
obtained at various q. The appearance of 0(7), as shown in Fig. |3l suggests the 
power law of the form ~ a(l — 7/70)^ -f b. Constants a and b can be found 
from the conditions that: (i) at 7 = 0 we have the stationary state with = 1 
for all individuals, hence 0(0) = 1; (ii) at 7 = 7 c the value of 0(7c) is given by 
Eq. m- Therefore, the only fitting parameter left is the exponent S. The fitting 
expression is 


0(7) = 


1-27c 
1 -7c 



( 10 ) 


and the results of fitting at selected q are shown in Fig. [3] via dashed lines. 

Therefore, the simulations of CA SIS model performed for various neigh¬ 
bourhood size q allow us to extend the zero dimensional SIS model to the case 
oi q < N. The solution for the fraction of infected individuals in the stationary 
state is given by (jH]). Here the analytic form (fTOll for 0(7) is found by fitting the 
simulation data obtained for CA SIS model. The dependence for the critical 
curing rate 7c and the exponent 6 on q can also be modelled from the simulation 
data, as shown in Fig. |4l Here the following power laws are used: 


Jc = l-Aq-\ 6 = Bq-\ 
where A ss 0.27, u ss 0.55, B Ri 0.56, and v ~ 0.25. 


( 11 ) 
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Figure 4: Dependencies of the critical curing rate 7 c and fitting exponent S in Eq. eg vs q 
for the extended zero-dimensional SIS model with fixed q. 


4. Cellular automaton SIS model with random neighbourhood size 

It has been pointed out that characterisation of each individual by its own 
intrinsic infectivity moves the description of the disease spread towards more 
realistic behaviour [l^. To do so we consider the CA SIS model, where the 
neighbourhood size qk is random for each fcth individual at each time step. The 
evolution of the system of N such individuals is run according to the algorithm 

• select random individual k 

• if Sfe = 0 do nothing 

• if Sfc = 1 then set Sfe = 0 with the probability 7, choose random neigh- 
bourood size qk from uniform distribution in [0; gmax] and set Sk' = 1 for 
one of its neighbours with the probability 1 — 7 

• perform N attempts described above to complete a single time-step 

This algorithm is interpreted as the one describing the set of individuals with 
random implicit mobility. The random neighbourhood size qk is globally bounded 
by 9max- As far as the critical curing rate 7c decays with the decrease of q (see, 
Fig.U]), its value for the model with random qk < gmax will be lower than that 
for the model with fixed neighbourhood size equal to (/max- Both models can be 
matched to have the same value of 7c, in which case one requires 

^ rqmakx. 

- / lc{q')dq' = -fdq) 

Q'max Jo 

where the l.h.s. represents the average value for 7c for the model with random 
neighbourhood size bounded by gmax; whereas the r.h.s. is the same value of 7c 
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Figure 5: The same as in Fig.[2]but for the CA SIS model with random neighbourhood size 
at various gmax- Values for critical curing rate 7 ^ are: 0.376 (gmax = 20), 0.426 (gmax = 36), 
0.462 (gmax = 88 ) and 0.491 (gmax = 420). 


obtained for the model with fixed neighbourhood size at particular q. Using the 
analytic fit for ■jdQ) dill) and performing integration we obtain simple relation 
between (/max and matching value of q 


^max Q 



4.279. 


( 12 ) 


The simulations performed at the values close to those suggested by the 
relation (HU): 9 max = 20, 36, 88 and 420 match the respective values for jc 
obtained for 9 = 4, 8 , 20 and 96 reasonably well, as shown in Fig. [5l 

One may conclude that in terms of the critical properties of the stationary 
state /(oo), the CA SIS model with random neighbourhood size can be matched 
well by its counterpart with appropriately chosen fixed value of 9 . 9 max and 9 
are related by a simple expression (IT^ obtained from the fit (ITTI) for 7 c vs q. In 
the next section we will focus the differences between these two models in terms 
of the spatial patterns and dynamics of approaching the stationary state. 


5. Spatial patterns and evolution of the system towards the stationary 

state 

One of principal advantages of using the CA model defined on a graph is the 
possibility to analyse spatial patterns of clustered infected individuals, as well 
as to trace their time evolution [HE El. In this study we analyse clustering 
effects or infected individuals both in the stationary state and when the system 
approaches it. To define clusters we introduce the cluster neighbourhood size, 
qc, which, in general, does not coincide with the infecting neighbourhood 9 or 
9 fe. For the case of considered here, the value 9 c = 4 is used (the nearest 














neighbours on the square lattice). Thereafter, the clusters are enumerated by 
index m and the cluster label Ck is assigned to each fcth individual having the 
meaning of the host cluster for this individual. The following algorithm for the 
clusters recognition is used: 

1 . set {cfe := 0}, m = 0; 

2 . pick fcth individual randomly with = 1 and = 0; 

3. set m := m + I, Ck := m, declare fc as a newcomer; 

4. loop over the cluster neighbourhoods {1} of all newcomers; 

5. if s; = 1 then set c; := m, declare I as a newcomer; 

6 . repeat steps 4-5 until no newcomers; 

7. go to 2. 

Alternatively, the Hoshen and Koppelmann [ 2 ^ algorithm can be used. 

At given curing rate 7 and time instance t we analyse the number of clusters 
^cil, t) and a set of their sizes {5^(7, t)}. Cluster size is given by the number of 
infected individuals it contains. Reduced number of clusters, as well as average 
and maximum cluster sizes at given t are given as: 


nc{'y,t) = 


acin,t) = 

^(*S'm(7i t))m, 

= 

^ma.x{Sm{7:t)} 


Their time averages performed in the stationary state are denoted as n-c {"/)■, 
acil) and me ( 7 ), respectively. 

The behaviour of these averages is shown in Fig.jSlfor both cases of CA SIS 
model with fixed q [frame (a)] and CA SIS model with random q [frame (b)]. 
The set of values q = 4,8,20,96 used for the model with fixed q matches the 
respective set of values Qmax = 20,36,88,400 for the model with random qt in 
terms of close respective values of 7 c, as discussed in Sec. IH Eq. (fT^ . For both 
types of models one can distinguish three following regimes: 

(a) single large cluster (7 < 7 +): nc{^) 0, 00 ( 7 ) > 0, 0.9 < rndj) < 1; 

(b) large -I- small clusters ( 7 + < 7 < 7 *): 0 < nc{l) < 2, 00 ( 7 ) 0, 0.25 < 

mci'f) < 0.9; 

(c) small clusters only ( 7 * < 7 < 7 c): neij) peaks, acil) ~ rnci'j) ~ 0 . 

The value of 7 ^ is defined approximately, whereas 7 * is the inflection point for 
the mcid- The regimes (a)-(c) are illustrated in the series of snapshots shown 
in Fig. [3 that are obtained for the smaller system of 60 x 60 individuals. 

The presence of these regimes indicates that besides the phase transition 
that occurs at 7 —^ 7 c, where I{oo) vanishes, there is another transition which 
takes place earlier, at 7 —>■ 7 *, where the size of the largest cluster of infected 
individuals decreases down to zero. At 7 * < 7 < 7 c, only small disconnected 
clusters are to be found [regime (c) above]. The position of this transition, given 
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Figure 6 : Reduced average cluster size ac( 7 ), maximum cluster size mci'y) and reduced 
number of clusters nci'y) vs curing rate 7 in the stationary state, (a) CA SIS model with fixed 
neighbourhood size q indicated in the figure; (b) CA SIS model with random neighbourhood 
size at various Qmax indicated in the figure. Characteristic values of 7 that separate various 
clustering regimes are denoted as 7 "*", 7 * and 7 c (critical curing rate). 



Figure 7: Snapshots with colour-coded clusters of infected individuals for CA SIS model 
with fixed neighbourhood size g = 4. Susceptible individuals are not shown, (a)-(c) illustrate 
respective clustering regimes described in the text. 
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Figure 8 : Time evolution of the reduced number of clusters nci'y^t) and the fraction of 
infected individuals for the CA SIS model with fixed neighbourhood size q indicated 

in the figure, (a) simulations at respective values of 7 that lead to /(oo) = 0.9 for each g; (b) 
the same for /(oo) = 0.7; (c) the same for I{oo) = 0.5. 


by 7*, is found to be weakly dependent on q for the CA SIS model with fixed q 
[Fig. [6] (a)]. However, despite matching both models via their respective values 
of 7c, their respective values of 7* do not match. This is especially evident 
for the cases of gmax = 20 vs q = 4 and gmax = 36 vs q = 8. In the model 
with random implicit mobility (^max = 20 and 36) the maximum cluster size 
vanishes at lower curing rate 7 when compared to the model with fixed implicit 
mobility (q = 4 and 8). Therefore, the randomness in the neighbourhood size q^ 
promotes stronger splitting of infected individuals into separate clusters. The 
higher maximum values of nc(n) occuring within the interval 7* < 7 < 7 c for 
the cases gmax = 20 and 36 as compared to the respective cases q = 4 and 8 also 
supports this conclusion. 

We will consider now the evolution of relevant properties towards the sta¬ 
tionary state in both models. To this end the series of runs are performed each 
started from the initial configuration containing a single infected individual. 
These series are split into three sets: (a) very low curing rate; (b) low curing 
rate; and (c) moderate curing rate. Each set contains a number of runs per¬ 
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Figure 9: The same as in Fig. 0 except for the CA SIS model with random neighbourhood 
size Qk < Qmax indicated in the figure. 


12 


























formed for both models at various values of q and (?max, respectively. The curing 
rate 7 is chosen individually for each run from the condition that the stationary 
state with /(oo) = 0.9, 0.7, and 0.5 is reached when the run belongs to (a), 
(b), and (c) set, respectively. The evolution of the instant fraction of infected 
individuals 7(7, t) and reduced number of clusters nc{'y,t) are plotted vs time 
t for both models in Figs. H] and IHl respectively. The solution ([S]) for 
obtained for the zero dimensional SIS model is also shown and is marked via 
q = N. 


Within each set, (a), (b), or (c), for the model with fixed neighbourhood 
size q, the increase of q leads to the reduction of a timescale for the system 
to reach the stationary state, monitored via behaviour of I{"f,t) (see. Fig. |S]). 
Such effect is to be expected, as far as the increase of q is interpreted in terms 
of enhancement of the implicit mobility, which promotes faster spatial spread of 
the infection. The fastest infecting dynamics is, obviously, achieved within the 
zero dimensional SIS model corresponding to the ideally mixed case SS110- 
Upon introducing random neighbourhood size qk, the infecting dynamics also 
speeds-up when compared to the results of respective matching runs performed 
at fixed q (compare Figs. [ 8 ] and [9l). The speed-up ranges from about 1.5 to 2 
times depending on the particular set, (a), (b), or (c). 

The evolution of spatial patterning also demonstrates marked dependence 
on the model type and the parameters of the run. In particular, for the case of 
fixed q and moderate curing rate [see, Fig.|Sl (c)], the reduced number of clusters 
n-ci'yjt) approaches its stationary value monotonically from below. However, 
the behaviour is different at very low and low curing rates [frames (a) and (b)]. 
In these cases, at larger neighbourhood sizes, q > 8, there is a pronounced 
maximum in the shape of nc{'y,t) curve, approximately a mid-way to the sta¬ 
tionary state. With the increase of q its height increases. The explanation for 
the presence of this maximum is the follows. At the early stage the system is 
mostly uninfected and for larger infection range the newly infected individuals 
will be formed at some distance from existing infected ones thus increasing the 
number of separate clusters nci'y,t). However, at a mid-way to the stationary 
state about half of the system is already infected and newly infected individuals 
begin to link the existing clusters, thus reducing the number nci'y,t). Simi¬ 
lar behaviour is also seen within the model with random qk (shown in Fig. |9]), 
where even higher maxima are observed when compared to matching cases of 
the model with fixed q. Therefore, the introduction of the randomness in the 
neighbourhood size enhances the hill-like shape in the evolution of the number 
of clusters nc{'y,t). At very low curing rate [frame (a)] the peaks are observed 
for all gmax being considered. 

Therefore, the dynamics of spatial patterning for infected individuals is found 
to depend strongly on the model type and the parameters of the run. When 
the initial state comprises almost healthy system (save for one individual), then 
for the contact-process-like infection spread, q = 4, the grows of a number of 
infected clusters towards its stationary value is always monotonic. With the 
increase of the neighbourhood size, g > 4, this behaviour changes and the 
number of clusters shows a pronounced maximum about a mid-way towards its 
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stationary value. This type of behaviour is getting more pronounced for the 
model with random neighbourhood size. 

6. Conclusions 

In the present work, we study the SIS model on geometric graph obtained 
from Z^. To each vertex the individual is attached with two possible states: 
susceptible or infected. Two options are considered for the neighbourhood size 
q, namely: the fixed q; and that taken at random at each step and for each 
individual qk globally bounded by gmax- The neighbourhood size is interpreted 
as the level of individuals’ implicit mobility. The evolution of the system of N 
individuals is run according to the asynchronous cellular automata algorithm. 

The stationary states are studied first, where we concentrated on the phase 
transition like behaviour for the ratio of infected individuals I{oo) which van¬ 
ishes at the critical curing rate 7 = 7 c. For the model with fixed g, with the 
increase of g, ■jc increases towards the value 0.5 characteristic for the zero di¬ 
mensional SIS model, which corresponds to the case of g = TV. The simulation 
data obtained at a range of g is used then to extend the solution for the zero 
dimensional model to cover the case q < N. Simple power-law model form is 
found for 7 c as the function of g. Close match is found for the phase transition 
like behaviours of the models with fixed g and random g^ bounded by gmax when 
g and gmax obey certain relation. 

In the stationary state we found three distinct regimes of spatial patterning 
of infected individuals depending on the curing rate 7 : (a) single cluster; (b) 
large-l-small clusters; and (c) small clusters only. The phase transition like 
change between (b) and (c) occurs at 7 *, the inflection point for the size of the 
largest cluster. The value of 7 * is found to be weakly dependent on g for the 
fixed neighbourhood size model but markedly dependent on gmax for the model 
with random neighbourhood size. It is found that randomness in g^ promotes 
splitting of the largest cluster at lower 7 as compared to the case of matched 
case with fixed g. 

Both the increase of the neighbourhood size and its randomness reduce the 
time needed by the system to reach the stationary state when simulation is 
started from the almost healthy state. This is understood in terms of the rela¬ 
tion between the neighbourhood size and the individuals’ implicit mobility. The 
type of evolution of spatial patterning also demonstrates marked dependence on 
the model type and the parameters of the run. For the contact-process-like in¬ 
fection spread, g = 4, the grows of a number of infected clusters towards its 
stationary value is always monotonic. With the increase of the neighbourhood 
size, g > 4, this behaviour changes and the number of clusters shows a pro¬ 
nounced maximum about a mid-way towards its stationary value. This type of 
evolution, which involves the initial built of a huge number of separated clusters 
and their following merge into larger ones, is getting more pronounced for the 
model with random neighbourhood size. 

This study, performed for a relatively simple SIS model can be extended 
towards more complex models involving latency, immunity, other types of graphs 
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and memory effects. 
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